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Abstract 

We consider the small mass asymptotic (Smoluchowski-Kramers approximation) 
for the Langevin equation with a variable friction coefficient. The friction coefficient 
is assumed to be vanishing within certain region. We introduce a regularization 
for this problem and study the limiting motion for the 1-dimensional case and a 
multidimensional model problem. The limiting motion is a Markov process on a 
projected space. We specify the generator and boundary condition of this limiting 
Markov process and prove the convergence. 
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1 Introduction 

The Langevin equation 

1*$ = &(<#) - \q? + a(q?)W t , q% = q € R n , q» = p € R n , (1.1) 

describes the motion of a particle of mass \i in a force field b(q), q £ M. n , subjected 
to random fluctuations and to a friction proportional to the velocity. Here Wt is the 
standard Wiener process in M n , A > is the friction coefficient. The vector field b(q) 
and the matrix function a(q) are assumed to be continuously differentiable and bounded 
together with their first derivatives. The matrix a(q) = ((Hj(q)) = cr(q)a*(q) is assumed 
to be non-degenerate. 

It is assumed usually that the friction coefficient A is a positive constant. Under 
this assumption, one can prove that q^ converges in probability as \x 1 uniformly on 
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each finite time interval [0, T] to an ra-dimensional diffusion process q t : for any k, T > 
and any p£ = p G R n , q% = q G W 1 fixed, 

limP I max \q+ — q+Ld > K ) = . 
/40 \0<kt' 1 UK J 

Here q t is the solution of equation 

Qt = \b{Qt) + \°{<lt)Wt , q = q% = 1 € M n . (1.2) 

The stochastic term in (1.2) should be understood in the ltd sense. 

The approximation of q^ by q t for < fi « 1 is called the Smoluchowski-Kramers 
approximation. This is the main justification for replacement of the second order equa- 
tion (1.1) by the first order equation (1.2). The price for such a simplification, in 
particular, consists of certain non-universality of equation (1.2): The white noise in 
(1.1) is an idealization of a more regular stochastic process W t with correlation radius 
5 « 1 converging to Wt as 6 J, 0. Let q^' S be the solution of equation (1.1) with W t 
replaced by W t . Then limit of q^' S as /x, 5 1 depends on the relation between fi and 
8. Say, if first 5 i and then fj, 4 0, the stochastic integral in (1.2) should be understood 
in the Ito sense; if first fi J, and then <5 4- 0, ,<5 converges to the solution of (1.2) with 
stochastic integral in the Stratonovich sense. (See, for instance, [5].) 

We considered in [6] the case of a variable friction coefficient A = A(q). We assumed 
in that work that X(q) is smooth and < Ao < A(q) < A < oo. It turns out that in this 
case the solution q^ 1 of (1.1) does not converge, in general, to the solution of (1.2) with 
A = A(q), so that the Smoluchowski-Kramers approximation should be modified. In 
order to do this, we considered in [6] equation (1.1) with W t replaced by W t described 
above: 

= b(q^) - A(q^)q^ + a(qf ^ , qf = q , qf = p . (1.3) 

It was proved in [6] that after such a regularization, the solution of (1.3) has a limit 
qf as fx | 0, and q\ is the unique solution of the equation obtained from (1.3) as fi = 0: 



4 ' = ^f) ^,( ' ^ ' ) + ^f) ^r< ^'•' !S = , • (1 ' 4) 



Now we can take 5 J, in (1.4). As the result we get the equation 

1 1 

= \{qtj b<yq ^ + \(q) a ^ ° Wt ' q ° = q ' ^' 5 ^ 

where the stochastic term should be understood in the Stratonovich sense. We have, for 
any 5, k, T > fixed and any p^ ,S = p fixed, that 
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and we have 



lim E max \qf — flf+Ld = . 
5^0 te[0,T] 



So the regularization leads to a modified Smoluchowski-Kramers equation (1.5). 

In this paper we study a further generalization of the problem considered in [6]. 
Keeping the assumptions on uniform boundedness and smoothness of A(«), we drop the 
assumption that < Ao < X(q) and instead assume that X(q) = for q £ [G] C MJ 1 and 
X(q) > for q £ Here G is a domain in M. n and [G] its closure in the standard 

Euclidean metric. For simplicity of presentation we assume in the rest of this paper that 
<j(«) is the identity matrix. (In Section 3 we further assume that &(•) = 0.) In order to 
use the results of j6] we introduce a further regularization of problem (1.5). We consider 
the problem 

and we study the limit of q\ as e I 0. This limiting process can be regarded as a limiting 
process of the system 

^ = b(q^) - (X(q^ £ ) + e)qf>* + W & t , qf' e = q , qf' e = p (1.7) 

as first n \. then 8 \, and then e I 0. 

System (1.6), in Ito's form, can be written as follows: 

However, as will be shown later, for non-compact region [G], it is sometimes more 
convenient to consider the projection of the above system onto another space X. (In 
particular, in Section 3 the space X is a cylinder X = S 1 x [a — 1, 6+ 1] for a < 0, 6 > 0.) 
Let us work with system (1.8) on X and compact region [G\. It turns out that, in the 
limit, to get a Markov process with continuous trajectories, one has to glue all the points 
of [G] and form a projected space £. Let the projection map be it : X — > <t. We will 
prove, for the 1-dimensional case (Section 2) and a multidimensional model problem 
(Section 3), that the processes (f t = Tr(qrf) converge weakly as e i to a continuous 
strong Markov process q t on C We will characterize the generator of this Markov 
process and specify its boundary condition. In particular, we will show that as e > is 
very small, certain mixing within [G] is likely to happen for the process qf . This mixing 
is the key mechanism that leads to our special boundary condition. We expect that (see 
Section 4), within the region that the friction is vanishing, similar mixing phenomenon 
will happen for the general multidimensional case. 

It is worth mentioning here that some related problems arc considered in [12J, 
[13], [15] and [IB] . It is also interesting to note that the limiting process for our two 
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dimensional model problem (see Section 3) shares some common feature with the so 
called Walsh's Brownian motion (see, for example [1]). 

However, at this stage we are not able to prove, in the most general multidi- 
mensional case (except for the 2-d model problem in Section 3), the convergence of 
qf = 7r(qf) in (1.8) to some Markov process q t . We will formulate a conjecture about 
this in Section 4. 



2 One dimensional case 

Let us consider in this section the 1-dimensional case. Besides the usual as- 
sumptions made in Section 1 we suppose that our friction A(») satisfies X(q) > for 
q € (— oo, — 1) U (l,oo). Let X(q) = for q £ [—1,1]. Equation (1.8) now takes the 
following form: 

qt ~ ~ 2(X(q £ t ) + e) 3 + XiqlH-e Wt ' 9 »-* £R ' (2J) 

We suppose that qo € [a — 1, b + 1] for some a < < b. The process qf is supposed 
to be stopped once it hits q = a — 1 or q = b + 1. 

Our goal is to study the asymptotic behavior of (2.1) as e 4 0. To this end we shall 
write the process (2.1) as a strong Markov process subject to a generalized second order 
differential operator in the form D v sD u s (see [1], [2], [H]). We have 

u £ (q) = f\\{x) + e) exp (-2 J* b(y)(X(y) + e)dy\ dx , (2.2) 

v £ (q) = 2 f\x{x) + e) exp ^2 J b(y)(X(y) + e)dy^j dx . (2.3) 

For fixed e > 0, the functions u £ and v £ are strictly increasing functions in their 
arguments. As e 1 0, they will converge uniformly on finite intervals to the functions u 
and v defined by 

u(q) = [" X(x) exp (-2 f b(y)X(y)dy] dx , (2.4) 
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v(q) = 2 A(x)exp ^2 J b(y)X{y)dyj dx . (2.5) 

The functions u and v are strictly increasing outside the interval [—1,1] and have 
constant stretches on [—1,1]. 

Consider a projection map it: we let 7r([— 1, 1]) = and n(q) = g + lforg<— 1 
and ir(q) = q — 1 for q > 1. Consider the process q £ = n(qf). Process qf for fixed e > 0, 
in general, is noi a Markov process. 



4 



Let us define two functions u and v as follows: u(q) = u(q — 1) for q < and 
u(q) = u(q + 1) for q > and u(0) = u(l) = u(— 1) = 0; = v(q — 1) for g < 

and t>(g) = u(g + 1) for g > and v(0) = v(l) = v(—l) = 0. Here the functions u and 
v are denned in (2.4), (2.5). The functions u and v are continuous strictly increasing 
functions on [a, b]. 

Define a Markov process qt on [a, b] as follows. The generator A of q\ is j4 = D^D^- 
The domain of definition -D(^4) of operator ^4 consists of all functions / that are con- 
tinuous on [a, b], are twice continuously differentiable in q G [a, b]\{0}, with finite limit 

lim Af(q) (taken as the value of Af(0)) and finite one-sided limits lim f ' ' ~ 



■I -r. ^ 4o u(S) - u(0) 

Dif(0) = Dzf( ) = lim M ~ f J-V 
uJy ' uJy ' sxo u(0) - u(-5) 

(taken as the value of Af (a) and Af(b)) 



D V(°) = D uf(°) = ^ -Sl f ~t % Also we have I imj4 /(«) = l im ^/(5) = o 

54.0 U(U) — U{ — 0) q^-a q->b 



Lemma 2.1. There exists the Markov process qt on [a,b\. 



Proof. The existence of such a process could be checked similarly as in Section 
2] . For the sake of completeness and comparison with results in the next section we shall 
check it here. To this end we use an equivalent formulation of the Hille-Yosida theorem 
(see Section 2] also [17} Theorem 2]). We check three conditions. 

• The domain D(A) is dense in the space C([a, b}). This is because we can ap- 
proximate every continuous function / with one that is constant in a neighborhood of 0. 
After that in the interior part of the intervals [a, 0) and (0, b], at a positive distance from 
0, with a smooth function. The approximating smooth function satisfy our boundary 
conditions since Af(0) = D±f(0) = D~f(0) = 0. 

• The maximum principle: if / G D{A) and the function / reaches its maximum at 
a point xq £ [a, b], then Af(xo) < 0. If xq ^ we have f'(xo) = and f"(xo) < and 

n n ti \ /"(so) u"(x ) , 

DM{X0) = *(so)«'(*o) " v>( X0 )(u>(x 0) y f {X0) * ° • 

If the maximum is achieved at 0, we consider the expansion 



f(x) = /(0) + D z f (0)(u(x) - 5(0)) + (Af(0) + o(l)) / (v(y) - v(0))Mv) ■ 

Jo 

The last integral is 0{u(x)v{x)) as x 0. Since D~f(0) > and D±/(0) < 0, by 
our boundary conditions at we get D^fiO) = 0. This implies that Af{0) < 0. 

• Existence of solution / G D(A) of A/ - Af = F for all F G C([a,b]). On 
each of the intervals [a, 0) and (0, b] the general solution of equation Xf — D^D^f = F, 
F G C([a,b]) can be written as 
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f ± (q) = f ± (q) + G ± (g) . 

Here satisfy the equation Xf^-D^D^ = F, / + (0+) = (or /~(0-) = 0), 

D±f+(0) = (or D~f-(0) = 0) and G ± (q) satisfy the equation XG ± - D^D^ = 0, 
G+(0+) = k+ (or G-(O-) = fcf), D±G+(0) = k+ (or D~G~(0) = k 2 ). Here jfef and 
k 2 are constants. Our boundary condition gives kf = k7 and k 2 = k 2 . The boundary 
condition D^D^f + {a) = D^D^f~{b) = singles out a unique / € -D(-A). □ 

We have 

Theorem 2.1. As e | 0, for fixed T > 0, i/ie process qf converges weakly in the 
space C[ 0j T]([a 5 &]) to process qt- 

The proof of this Theorem is based on an application of the machinery developed 
in [51 Ch.8], [S] and [7]. We shall use the following lemma, which is the Lemma 3.1 of 
[HJ Ch.8, page 301]. We formulate it here in the terminology that meets our purpose. 

Lemma 2.2. Let M be a metric space; Y , a continuous mapping M i— > Y(M), 
Y{M) being a complete separable metric space. Let (Xf,P^.) be a family of Markov 
processes in M ; suppose that the process Y{Xf) has continuous trajectories. Let (yt, P y ) 
be a Markov process with continuous paths in Y{M) whose infinitesimal operator is A 
with domain of definition D(A). Let T > 0. Let us suppose that the space C[q^(Y(M)) 
of continuous functions on [0, T] with values in Y(M) is taken as the sample space, so 
that the distribution of the process in the space of continuous functions is simply P y . 
Let ^ be a subset of the space CiQ t00 \(Y(M)) such that for measures \xi on Y{M) 

the quality J F dn = J ' Fi n for all F £ * ,m P l,e, n = „, Let D h the suiset of 

D(A) such that for every and A > the equation Xf — Af = F has a solution 

fGD. 

Suppose that for every x £ M the family of distributions Q| ofY(X^) in the space 
CrQ t T](Y (M)) corresponding to the probabilities of P| is weakly pre-compact; and that 
for every compact K C Y{M), for every f € D and every X > 0, 

roo 

E% / e^[Xf{Y{Xt)) - Af{Y(Xf))]dt -> f(Y(x)) 
J o 

as e \. uniformly in x G y~ 1 (X). 

Then Q| converges weakly as e \. to the probability measure Pyw. 
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Proof of Theorem 2.1. Making use of Lemma 2.2, we take the metric space 
M = [a — 1,6 + 1] and the mapping Y = ir. The space Y(M) = ir([a - 1, 6 + 1]) = [a, 6]. 
We take the process qf as (X^,P £ X ). We take the process q\ as (yt,P y ). 

Let \T/ be the space of all continuous bounded functions in [a, b] which are once 
continuously differentiable inside [a, 0) and (0,6], with bounded derivatives. The space 
D C D(A) consists of those functions / G D(A) such that they are continuous and 
bounded in [a, b] and are three times continuously differentiable inside [a, 0) and (0,6], 
with bounded derivatives up to the third order. 

Pre-compactness of the family of distributions of the process {ql} e >o is checked in 
Lemma 2.4. What remains to do is to check that for every compact K C [a, 6], for every 
/ G D and every A > 0, 



as e \, uniformly in go € vr 1 (K). This is done in Lemma 2.5. This finishes the proof 
of Theorem 2.1. □ 

For positive 6 small enough, let G(5) = [a- 1, -1 - 6] U [1 + 6, 6 + 1]. Let < 5' < 5. 
Let C(6') = {—1 — 6', 1 + 5'}. We introduce a sequence of stopping times To < ao < t\ < 
01 < t 2 < 02 < ••• by 

r = , a n = min{t > r n ,qf G G(5)} , r n = min{t > a n _i : qf G C(5')} . 

This is well-defined up to some at (k>0) such that 

~Pq% k (Qt+a h hits a — 1 or 6 + 1 before it hits — 1 — <5'orl + <5) = l. 

We will then define Tk+i = min{t > at '■ of = a — 1 or 6 + 1}. And we define 
Tk+i < o-fc+i = Tfc+i + 1 < T k+2 = T k +i + 2 < a k+2 = Tfc+i + 3 < ... and so on. 

We have lim r n = lim a n = oo. And we have obvious relations q% G C(5'), 

n— >oo n— >oo " 

q% n G C(S) for 1 < n < k (as long as k > 1, if k = the process may start from G(5) 
and goes directly to a — 1 or 6 + 1 without touching C(<5') and is stopped there, or it 
may start from (—1 — 5, 1 + 5), reaches {—1 — 5,1 + 5} first and then goes directly to 
a — 1 or 6 + 1 without touching C(5') and is stopped there). Also, for n > k + 1 we have 
q e Tn = q £ an = a — 1 or 6 + 1. If q^ = qo G G(<5), then we have do = and T\ is the first 
time at which the process qf reaches C(5') or {a — 1, 6 + 1}. 

Now we check weak pre-compactness of the family of distributions of the processes 
{qf}e>o- To this end we need the following lemma, which is Lemma 5.1 in [7J. We 
formulate it using our terminology. 



E, 




[A/(7r(gf )) - Af(7r(qf))]dt - f{v{q )) 
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Lemma 2.3. Let (f, ,S for every e > 0, 5 > 0, be a random element in C[ ,T]([ a ) b]) 
such that maxjgjr — q^' S \ < 5 on the whole probability space. If for every positive 5 the 

family of distributions of ql' S , e > 0, is tight, then the family of distributions of q% is 
pre- compact. 

Now we have 

Lemma 2.4. The family of distributions of {ql} £ >o is pre-compact. 

Proof. Let 5' = 8/2 so that we need only one parameter <5. Between the times Oj_i 
and Tj the process qf is either in [a,— 1 — 5/2) or in (l + <5/2, b], and for cjj_i < t < t' < T{ 
we have |g| — gjr, | = |gf — gf, |. Since we have 

qt ""-I [K q J)T-e-2(X(q!) + e^\ ds + J t HqjH-e dW - 
we can estimate 

E\q!-q!,\ 4 <K(5)\t-t'\ 2 . 
The constant K{8) is independent of e provided that e is small. Now we let 



- y o i GW2) (5.) [t^^ - 2(A(9f)+e)3 j d * + y ^(o^^. . 

From the above estimate we see that Z £ t )& for fixed <5 > is tight. The trajectories 
of these stochastic processes satisfy the Holder condition \Z £,S — Z £ t f\ < H £ ' S \t — i'l 1 / 5 
where H 6 ' 5 are random variables with E(i? £,<5 ) 4 bounded by the same K(8). 

For i > 1 if q\. G C(<5/2) and g £ . G C(<5) then between the times Tj and <7j (< T) 
the process gf travels a distance at least (5/2 and at least this distance in (7(5/2) on the 
same interval either [a, —1 — 5/2) or (1 + 5/2, b\. By our estimate on Holder continuity 

of Z £,S this implies that <jj — Tj > f ^ ^ , i > 1. If g £ . G {a — 1, 6 + 1} then by our 

definition of the stopping time <7j = r% + 1 we can choose 5 small enough such that the 
above inequality also holds. 

Now we shall define the process q^ as follows. 

• For <7j_i < t < Ti we take q^' S = qf. 

• For To < t < do we take gf' 5 = gL . This gives max |g^ ,<5 — gf | = max |gL — 

ro<t<cro ro<t<ao 

ql\<5. 

• If Tj < T < ai we take g^' 5 = g^ for n <t <T. This gives max |g^ ,<5 — gf | = 
max 125. — gf I < (5/2. 

r t <t<T 1 r ' 1 ~~ 
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• If Oi < T. In this case if q e Ti and q% % are within a distance < S from 0, we define 

Q T i+ ^ = 0, 

2 

~e,S f, 2(t-Tj) \ ^ Ti + Oi 

i t ={ 1 --^r)««*xii<t<-j-, 

Qt =-{l-^-^-) qa Jor^<t<* t . 
Since this is just a linear interpolation it is clear that in this case we have max \(f t ' S - 

Ti<t<Ui 

gf | < 25. Within this time interval Tj < t < t' < Oi, i > 1 we have 

|^ _ < 6 A t - t'\ < ^ \t - t'f/ 5 < 2 u ^H £ ' S \t - if/ 5 . 

Wi Ti| (min-K -n\y/s 

i>l Z 

Another possibility is that q%. = q%. = a — 1 or 6+1. In this case we define q^' S = qf 
for Ti < t < Oi. 

On the whole interval < t < t' < T we have \q/ - q^/\ < (2 11 / 5 + 2)H £ ' 5 \t' - 1] 1 / 5 
5 ^ 5 



for \t' — t\ < I - x . This means that for fixed 5 > we have the tightness of the 

family of distributions of q^' S in the space C[ y]([a, 6]). Since we checked maxj^'' 5 — 
qfl < 25, by using Lemma 2.3 with 25 instead of 5 we get the pre-compactness of the 
family of distributions of gf in C[ ,T]([ a > b]). □ 

The proof of the next Lemma 2.5 is based on Lemmas 2.6-2.10. Within the proof 
of this lemma and the auxiliary Lemmas 2.6-2.10, we will take e I 0, 5 = 5(e) I 0, 
5' = 5'(e) I in an asymptotic order such that < e « 5' « 5. Although not very 
precise, but for simplicity of presentation we will just refer this choice of order as first 
e I 0, then 5' I and then 510. It could be checked that such an order of taking limit 
does not alter the validity of the result. 

Throughout the rest of this section and next section when we use symbols U, V, 
Mi, Ci, Ai, etc., they are referring to some positive constants. We will not point out 
this explicitly unless some special properties of the implied constants are stressed. Also 
we sometimes use the same letter for constants in different estimates. 

Lemma 2.5. For every compact K C [a, b], for every f G D and every A > 0, 

e- Xt [\f(7r(q!)) - Af(ir(qf))]dt - /(vr(g )) -> 
as e 10 uniformly in qo £ -n~ l (K). 

Proof. The above expectation can be written as 
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E 



qo 



E 

n=0 



A<T n 



+ 



r r T "+ i 



n=0 L - y °"" 

oo 



E. 



qo 



where 



.71=0 



n=0 



(2.6) 



^i(?) = E 9 



/ e^[Xf(n( q f))-Af(n( q mdt + e-^f(n( q ^)) - /(*(*)) , (2.7) 
y 



/ Tl e-^[A/(7r(gD) - A/(7r(gf))](ft + e-^/Mfc)) 



-/(*(?))■ (2-8) 



7 £ ^ 



We used the strong Markov property of qf. Since for n > k + 1 we have ipKq^ 
ip^{q% ) = we can assume that the function ^| is taken at a point on G(5)\{a — 1, 6+1} 
and the expectation is determined by the values of the process qf in one of the intervals 
either (1 + 5', b + 1] or [a — 1, —1 — 5'). We will prove, in Lemma 2.6, that under our 
specified asymptotic order we can have IV^Oz)! — (^(^) ~~ ^)) 2 as £ i 0- 

We can assume that the function V>f is taken at a point in [—1 — S',1 + 5'] (in the 
case when n = and g§ G G(<5), we also have tpf (<7o) = 0). We can write 



= (E <? /(vr(^ ))-/(^))) -E,(l-e- A -)/(vr(^ ))+E 9 y o e -*[A/(7r(q?)) - A f (*{<£)))]& 
= (iy(q) + (Iinq) + (IIiy(q). 

(2.9) 

We are going to prove, in Lemma 2.8, that for q £ [—1 — 5', 1 + 5'], for a function 
/eflwe can have the estimate |(i") e (g)| < M\(u{5) — u{—5)) 2 . 

In Lemma 2.9 we will show that E g cro < Mi(u(<5) — u(— 6))(v(6) — v(—6)) and 
E,(l - e- Xa ") < Mi(iZ(<5) - 5(-($))(v(($) - 5(-<5)) so that |(//) £ (g)| + |(///) £ (g)| < 
Mi(5(<5) - 5(-<5))(5(<5) - «(-£)) for g € [-1 - 5', 1 + 5']. 

These estimates show that 

< (u(5) - u(-5)) 2 + Mi (5(5) - 5(-<S))(u(<5) - «(-<*)) 
for all g G [-1 - 5', 1 + <5']. 



10 



As we only consider the arguments q% n of ipf in (2.6) being in [—1 — 6', 1 + 5'] starting 
with n = 1 (otherwise ip[ = 0), we have, by strong Markov property of qf, that 



e, J>~ Arn V4« 

n=l 

< (5(5) -u(-6)) 2 + M 1 (u(6) -u(-6))(v(6) -v(-6))J2^~ XTn 



n=l 

oo 



n-1 



< (5(5) - 5(-5)) 2 + M l (5(5) - 5(-5))(5(5) - v{-6)) V sup E ge - Ari 

^ V^G(«5) / 
We will show, in Lemma 2.10, that E 9 e- AT1 < 1 - M 2 5(5) A (-5(-5)) for all 



q € G(5). Since as 5 I we have < M 2 < 



5(5) 



-5(-5) 



< M3 < 00, we have 



n=l 



< ((5(5) - 5(-5)) 2 + M x (5(5) - 5(-5))(5(5) - 5(-5)))- 



M 2 (5(5)) A(-5(-5)) 

as 5 4- 0. For n = the expectation E go ^f (gg) is small as e is small. 
For the second term in (2.6) we can estimate 



n=0 
<(1 + 



M 4 



< ^E g e- A -|^(g)| < ^E 9 e- A -|VK<?)I 

n=0 n=0 

)(5(5)-5(-5)) 2 



(5(5)) A (-5(-5)) 
which converges to as e 4- 0. This proves this lemma. □ 



Lemma 2.6. We have, for q € G(5), as e is small, that IV'iWI — ~~ u (~ $)) 2 ■ 

Proof. For the initial point q € G(5) and the time interval < t < t\ the trajectory 

of q\ is traveling in one of the intervals either [1 + 5', 1 + b] or [a — 1, —1 — 5']. Without 

loss of generality let us assume that q € [1 + 5, 1 + b] and we are traveling in the interval 

[1 + 5', 1 + 6]. Let q = ir(q). Let B(q) = b{q + 1) and A(q) = X(q + 1). Let us extend the 

function A(«) to the whole line R. The extended function A(«) is smooth, bounded, with 

uniformly bounded derivatives and such that A(x) > min A(g), A(x) = A(l + x) 

q e[i+S',i+b] 

for x e [6',b\. 

Let the process q t be subject to the stochastic differential equation 



Qt 



A'(f) 



+ 



A(f t ) + e 2(A(f)+e)3 l(q t ) + e 



W t , q = q ,0 <t < 00 
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We introduce a stochastic process q t , q = q with generator A, subject to the 
stochastic differential equation 



Qt 



B{q t ) A'(q t 



+ 



1 



■W t , q = q ,0 < t < oo . 



A(q t ) 2A*(q t ) A(q t ) 
Notice that the modified generator A agrees with A before the process gf reaches 



q^ ri . And before the time t\ the process q t agrees with the process gf. Therefore we 
have, 



Mi) = E ? 



n 



e- Xt [Xf(f t )-Af(q t )]dt-e-^f(q 



-An 



It is clear by Ito's formula that we have (also see, [101 Section 2]), for the stopping 
time tl, 



e- M [\f(q t )-Af(q t )]dt-e- X ^f(q 



Notice that the function / G D C D(A) is three times continuously differentiable 
in [<$', b]. This gives the estimate that for some positive U, V > and T = T{e) we have 



1^1(9)1 



% / e- xt [X(f(f t ) - f(q t )) - (Af(f t ) - Af(q t ))]dt - e^UO " /(5n)) 



T(e) 



Ae- A 'dt(Lip(/)).|5;-^| + 



jf e- At dt (Lip(A/)) • |gj - + (Lip(/)) • |r n - gjl (ri < T(e)) j + 
VP(n > T(e)) 



< C/ ( 0< max £) - q t \ J + VP(ti > T(e)) 

< U^jdsx (Eq\f t -l\ 2 ) 1/2 + VP(n > T(e)) . 



By the integral form of the stochastic differential equations of the processes q t and 
q t we have 
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\Qt~Qt\ 2 
< C 



f 




Jo 


[( 



B(f a ) 



A'(C) 



'B(Z) fl'(f s ) 



A(q s ) + e 2(A(q s )+e)\ 
(B{f s ) A'g) \ _ f B@ 8 ) 



Xg) 2(A(^)) 3 / 



r/,s 



+ 



A'(g s 



+ 



A(g) + e A{f s ) 



dW, 



+ 



A(S) A(gJ 



dW, 



Let a(A) be the Lipschitz constant of — (x > A), /3(A) that of — ^ (x > A), 7(5') 
^ x 2ar 

that of 2^ - -AM- (g > <5'), n(8') that of ^J— (g > 8'). Let m(<5') = min A(x). 
A(q) 2A(g) 3 A(g) xe[8',b] 

We can estimate 



A'(?) 



A(g) + e 2(A(^)+e)3 
<4 1 (tV[a 2 M5')) + ^( m (i'))]) , 

(B(f s ) A'(f s ) \ 



( B{f s ) 



/ 



/ VA(C) 2(A(S)) 
A'(? s ) 



(is 



\A(t) 2(A(t)) 3 KHZ) 2(A(f.))S 



ds 



; (<5') / E g ~|^ - 

JO 



g s | ds , 



A(E)+e A(J) 



dWs 



< / eV(m(<5'))ds = e 2 ta 2 (m(<J')) , 
Jo 



A(C) A(gJ 



</. 2 (<5') / E^-?/ds . 
Jo 



We have, by using the above estimates, with a possible change of the constant C, 



that 



V«\Tt-Qt\ 2 < C (te 2 (t(a\m(5>)) + f3 2 (m(5'))) + a 2 (m(6>))) + (t 7 2 (8') + fi 2 (6')) J* E^ -fj 2 ds 
By Bellman-Gronwall inequality we have 

Vq\Tt ~Qt\ 2 < Cte 2 (t(a 2 (m(5')) + f3 2 (m(5'))) + a 2 (m(5'))) exp (C(t 7 2 (0 + ^ 2 (<5'))*) • 
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1 A3 A% 

As we can check that \a(m(S'))\ < — ttftv, P( m (^')) < — ittttj 7(^') < — jtttt an d 



to 2 ((5') ' ' " - m 4 (<5')' ,v ' ~ m\&) 

this crivps as A' is smn 



A 3 

|m(<5')I < — 0/ r <s ) this gives, as 8' is small, that 



1/2 



CT(e)e(a 2 (m(5')) + /3 2 (m(<*')) + ^0Ay/2 exp ( C (T(e) 7 2 (5') + Ai 2 (<S'))T(e)) 



<CT{e) . £ exp|cT 2 ( £ ) . 

mm A 4 (g) I mm A 8 (g) / 

g6[l+5',l+6] \ qe[l+S',l+b] ' ' 

Noticing that by strong Markov property P(ti > T(e)) < K exp(— pT(e)) for some 
p > 0, K > 0, we see that 



qe[l+5',l+b] \ qe[l+5',l+b] 



Lei us choose '/'( ;) - \ In In-. We will then have 



/ 1\V 2 £ / min C A« ( 9 ) / f 

|^(<Z)| < C In -J — X%) V n I J < ' e[1+5 '' 1+61 + V r exp(- P yinln-)) . 



?6[l+5',l+6] 

For fixed 5' > 0, one can choose e small enough such that 



g6[l+(5',l+6]U[-l+o,-l-5'] 

for some Uq > 0, p > and < k < 1. As we choose first e | and then <5' J, 0, this 
gives that as e is small we have IV'IWI ^ ("(^) ~~ u(—5)) 2 . □ 

Lemma 2.7. M^e have, as £,5,6' are small, for q G [—1 — 5', 1 + <5'] and C > 0, 

that 



5(<5) - 5(0) 



P,(7r(gL) = -*) 



5(g) - 5(0) + Cg 
- 5(5) - 5(-<5) ' 

5(<5') - 5(0) + 



u{5) - u(-5) 



< 



u(S) - u(-S) 
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Proof. Let q = n(q) G [— 8', 6']. We have, for bounded positive functions C\(8, e), 
C2(8, e) and positive constants C\, C2, C, that 



PMqL) = 8) - 



5(0) - u(-8) 



< 



< 



dao/ ' u{8)-u{-8) 
u £ {q)-u £ {-\-8) 5(0)-5(-<5) 



u e {l + 8)-u £ (-l-8) u(8)-u(-8) 
5(0) - u{-8) + 5(5) - 5(0) + d (8, e)e ' 5(0) - u{-8) 



u{8) -u{-8) + C 2 (8,e)e u{8) -u{-8) 

(5(5) - 5(0) + C l£ )(u(8) - u{-8)) + C 2 e(5(0) - u{-8)) 



u{8') - 5(0) + Ce 



(u(8) - u{-8)f 



u{8) - u{-8) ' 
The estimate of P g (7r(g^ ) = —8) is similar. □ 



Lemma 2.8. We have, as e are small, for q € [—1 — 8', 1 + 8'], that \(I) £ (q)\ < 
C{u{8)-u{-8)) 2 . 



Proof. We have, using Lemma 2.7, that 



\{If{Q)\ 

= |E,/(7r(^ ))-/(7r(g))| 

(f(8) - /(0))P ? (tt(^ o ) = 8)- (/(0) - /(-<5))P ? (7r(^ ) = -8) + (/(0) - /(tt(q)))| 



< 



- /(0)) 5(5)-5(-5) " (/(0) " f{ - 5)) u(8)-u(-8) 
u(8>) - 5(0) + Me 

Q 5(5)-5(-5) + Cs(u(<5 } " U{0)) 

(u(Q) - 5(-5))(5(5) - 5(0)) f f(8) - /(0) /(0) - /(-*) 

5(<5) - u(-6) \u(8) - 5(0) 5(0) - u{-8) 

C 4 5( ? / : 5( ? / + , M£ + C 5 (5(^)-5(0)) 



+ 



+ 



5(5) - 5(-<5) 



< C 3 (u(8)-u(-8)Y + C 4 - 



5(5) - 5(-<5) 



+ C 5 («(<5')-«(0)) • 



We have used our gluing condition D~f(0) = D~f(0). Now we choose first e | 
then 5' I 0, we get, as e is small, that \{I) e {q)\ < C(u(8) - u(-8)) 2 . □ 

Lemma 2.9. As e, 8, 8' are small, for q € [—1 — 8', 1 + 8'] we have, 



B q a < C(u{8)-u(-8)){v(8)-v(-8)) , E,(l- e - Aff °) < C{u{8)-u{-8))(v(8)-v(-8)) . 



15 



Proof. We apply the well known formula for the expected exit time (see, for 
example [LA Chapter VII, Theorem 3.6]) and we have 

rl+S 

E g£ x = / G £ (q,r)dv £ (r) , 
J-l-8 

where the Green function 

( (u £ (q) - u £ (-l - 5))(u £ {l + 8)- u £ (r)) 

u £ (l + 8) - u £ (-l - 8) 
(u £ (r) - u £ (-l - 8))(u £ (l + 8)- u £ {q)) 



G £ (q,r)={ 



u £ (l + 8) - u e (-l - 8) 
otherwise . 



for -1- 5 <q <r <1 + S , 
for -1-8 <r <q<l+8 , 



Therefore it is easy to estimate 

< (u £ (l + 8)- u £ {-l - 8))(v £ (l + 8)- v £ (-l - 8)) 

< (u(8) - u(-8) + C 6 e)(v(8) - v(-8) + C 7 e) 

< C(u(8) - u(-8))(v(S) - v(-S)) 

as desired. 

This helps us to find 



E q (l-e- Xa °) = AE g 



< XB q a < C(u(8) - u(-8))(v(8) - v(-8)) 



□ 



Lemma 2.10. For q £ G(8) and 8 sufficiently small, we have 
limlimE g e- A ^ < 1 - C(u(8)) A (-«(-<*)) . 

5' 10 eiO 



Proof. Without loss of generality let q G [1 + 5,1 + 6]. The expected value 
M £ (q) = E g e~ AT1 is the solution of the differential equation D v eD u eM £ {q) = XM £ (q), 
M £ (l + 8') = M £ (l + b) = 1. 

There exist two solutions fi{q), f^il) °f the equation D v D u f = Xf with / A (l) = 
+ b) = 1 and / X A (1 + b) = / 2 A (1) = 0. The derivatives D u f?(x), D u f£(x) are 
increasing functions, — oo < limD u (f^ + f£ )(q) < 0, < lim D u (f£ + f^il) < oo (see 

94-1 9tl+6 

ID, HE). 

We shall make use of Lemma 2.6. Since q € [1 + 5, 1 + 6] we see that do = 0. Lemma 
2.6 tells us that, for k = 1,2, we have 
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lim 

£4-0 



E„ 



I" e-^Xf^qf) - D v D u f£(q!)]dt + 
Jo 



< { u{5)-u{-5)f 



Taking into account the definitions of / A , f£ we see that the above inequality gives 



limE ? e-^V fe A (^) -/£(<?) 



< (u( S ) - u{-5)f 



Since = f£(l + 5') when q e Ti = 1 + 5' and = / fc A (l + 6) when = 1 + 6, 

we see that for some K > we have 



limE e-^ - (/2 (1 + b) - / 2 A (1 + 5'))^) + + *') - #(1 + 6))/ 2 A (g) 

^ 9 / A(l + (50/ A ( l + 6) _ / A (1 + 6)/ A (1 + ( 5/ ) 

(The expression 

(/ 2 A (1 + b) - / 2 A (1 + + + g) - A A (1 + b))/ 2 A (g) 

/ A (l + 5')fi(l + 6) - / A (l + 6)/ A (l + 6') 

is the solution of the equation \f(q) = D v D u f with f(l + <5') = f(l + b) = 1.) 
This gives 



< K(S(<5)-£(-<5)) 5 



limhmE ? (l - e "^) - [1 - (/ A (g) + / A (g))] 
<5'4-0 £j,o 

pA i fA\ 



< K{u(5) -u{-5)) 2 . 



Taking into account that -oo < limD u (/ 1 A + / 2 A )(q) < 0, < lim D u (f? + f£)(q) < 
oo we see from the above estimate that 

limlimE ? (l - e~ Ari ) > C(u(5)) 

S'lO eiO 

for q G [1 + 5, 1 + b] and (5 sufficiently small. The case of u(— S) is handled in a similar 
way. □ 



3 A two dimensional model problem 

In this section we discuss a two dimensional model problem. We work with a 
Smoluchowski-Kramers approximation in the plane M 2 . Let us suppose that the friction 
coefficient A(«) depends on the y variable only: X(x,y) = X(y). Suppose for y € [—1, 1] 
we have X(y) = 0. For y £ [— 1, 1] we have X(y) > 0. For simplicity of presentation we 
also assume that the drift is zero: = 0. All the other assumptions about A(») are 
the same as was made in Section 1. 
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In addition, we assume that for e > 0, 

~1 1 /•!+£ 



(In the case that both integrals converge the proof of Lemma 3.1 repeat that in the case 
of both integrals divergent but we do not know anything about the case of one integral 
convergent and the other divergent.) 

As we already introduced in equation (1.8) of Section 1, we are actually considering 
the stochastic differential equation for the position of the particle q\ € IR 2 as follows: 

By taking into account our assumption on the friction coefficient A we can write 
the above equation in coordinate form. Let q\ = (xf,yf). Let Wt = (W^,W 2 ). We 
have 

'* ! -A5S+i* ?ia «- a * 6R ' (32) 

Let a < < b be given. Throughout this section we will assume that our process 
q\ is stopped once it exits from the domain {(x,y) 6 M. 2 : a — 1 < y < 6+1}. We 
therefore suppose that yo £ [a — 1, b + 1]. 

Note that, similarly as in Section 2, the process yf is a strong Markov process 
subject to a generalized second order differential operator in the form D v e^D u e^ where 



u £ (y) 

Let 



fW) + e)ds , v £ {y) = 2 F {\{s) + e)ds . (3.3) 
Jo Jo 



u{y) = [ V X(s)ds , v(y) = 2 f X(s)ds . (3.4) 
Jo Jo 

We have the obvious relation u e (y) = u(y) + ey and v £ (y) = v(y) + ley. 

Let us identify points in the x direction x ~ x + 2tt. Therefore we get a process on 
the cylinder S* 1 x [a — 1, b + 1], stopped once it hits the boundary {y = a — 1 or b + 1}. 
Let 

J 0f = xf mod 2tt , 
1 Vt = Vt ■ 

In the rest of this section we refer to the process q\ as the one on a cylinder: 
Qt = iPtiVt) ls on the cylinder S 1 x [a — 1,6 + 1]. When we speak about the process 
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q\ on the domain {(x, y) G IR 2 : a — 1 < y < b + 1} C M 2 we will instead refer to the 
coordinate representation (xf,yf). 

Let (£ be the product S 1 x [a, b] with all points S 1 x {0} identified, forming the 
point o. A generic point on <£ will be denoted q = (9,y) where 9 £ S 1 and y G [a, 6]. All 
points (6,0) correspond to o. 

Let us consider the following projection map 7r : S 1 x [a — 1, b + 1] — > <£. We let 

' , for Ky<6+1 ; 

7r(^,y) = < ! + fora-l<y<-l; (3.5) 

, for — 1 < y < 1 . 

Let 7r(qf) = qf = (#f,yf). We see that yf = -ir(yf) where tt is the projection map 
introduced in Section 2. 

Let, as in Section 2, u(y) = u(y — 1) for y < and u(y) = u(y + 1) for y > and 
5(0) = u(l) = u(— 1); iJ(y) = v(y — 1) for y < and ?;(y) = -u(y + 1) for y > and 
v(0) = v(l) = v{— 1). The functions S(y) and u(y) are continuous strictly increasing 
functions on [a, b]. Let X(y) = X(y — 1) for y < and A(y) = A(y + 1) for y > and 
A(0) = 0. 

Let A be the operator given, for y ^ 0, by the formula 

Af(0, y) = D m D m f + ~^fiif • (3-6) 

Let D(A) be the subset of the space C(£) consisting of functions /(q) for which 
Af(9, y) is defined and continuous for y ^ 0, the derivatives in it being continuous; such 
that finite limits 

, JjP _ B ■ „J™ 0+ Al.S./C', S) ■ (3-7) 



exist; 



exists and does not depend on 9; 



lim Af(9',y) (3.8) 
fl'->-e,i/-5>o 



lim A/(0',y) = lim A/(0',y) = ; (3.9) 

6>'->-6>,j/-s>a 6'^6,y-^>b 



and 

r 2n i-2-K 



L ^- D ^ !{e '^ de = L r 32*» D *a f(!f < w » ■ (3 - 10) 

It is worth mentioning here that the above condition (3.10) in the definition of D(A) 
can be replaced by the condition that lim Dfi(y)f(0', y) and lim D^y^f(6', y) 
not depending on 9 and coinciding. In this case the proof of Lemma 3.1 remains the 
same. 
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Let us define, for / G D(A), Af(9, a) and Af(0, b) as the limits (3.9) and Af(o) as 
the limit (3.8). The operator A defined on D(A) is a linear operator D(A) C(€). 



Lemma 3.1. The closure A\ D ^ of the operator A\ D ^ exists and is the infinites- 
imal operator of a Markov semigroup on C(C). 

(The corresponding Markov process q t stops after reaching the boundary of (£ 
(y = a or b).) 

Proof. We use the Hille-Yosida theorem and we check the following: 

• The domain D(A) is dense in C(£). 

This is because we can approximate every function g in C(C) by a function / which 
is smooth, close to g outside a neighborhood of o and is equal to g(o) in the neighborhood 
of 0. This function / satisfies our restrictions on D{A) and can approximate the function 
g with respect to the norm of C(£) as we choose the neighborhood of o small enough. 

• The operator ^4|d(A) satisfies the maximum principle: for / G D(A), if this 

function reaches its maximum value at a point q G £ we have Af(q) < 0. 

Indeed, for q = (9, a) or (9,b), we have Af(q) = 0. If q = (9,y), y ^ the first 

d 2 

partial derivatives at q are equal to and 7^2 /(^y) < 0) ^v{y)^u{y) — 0- Finally, 
if q = we have the left-hand derivative D~,~.f(6,Q) > 0, the right-hand derivative 
D~^~.f(9, 0) < and by (3.10) both these derivatives are equal to 0. It follows then that 
the limit as y — > of the second y-derivative is non-positive for all 9 G S 1 . Since the 
integral over S* 1 of the second 9 derivative is equal to for all y ^ 0, taking into account 
that Af{o) is equal to the limit (3.8), we have that Af(o) < 0. 

It follows from the maximum principle that for A > the operator XI — A\ D ^ 
does not send to zero any function that is not equal 0, and this linear operator has an 
inverse (that is not defined on the whole C ((£)), with ||(A7 — A\ D ^)^ 1 \\ < A -1 . Every 
bounded linear operator does have a closure (which is just its extension by continuity), 
and with it the operators XI — A\ D r A \ and ^Iiva) also have closures. 

• Finally, to check that we can apply Hille-Yosida theorem to the closure ^4|d(A) 
we have only to check that the bounded operator (XI — ^Iz^/l))" 1 is defined on a dense 
set. That is, for a dense subset of F G C(£) there exists a solution / G D(A) of the 
equation 

Xf-Af = F. (3.11) 

Let us take F(9,y) = e md G(y), defining F(o) as its limit as y — > 0. Of course for 
n/Owe have to have lim G(y) (which limit we'll take as the value G(0)) equal to 0. 

We shall look for the solution / G D(A) of the equation (3.11) in the form f(9,y) = 
^ in9 9{y) (again, for n ^ it should be g(0) = lim g(y) = 0). 
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The differential equation for g(y) following from (3.11) is the ordinary differential 
equation 

rr „ „ „ . _ (3.12) 



A 2 (y) 



and it should be solved with the boundary conditions 



n 



A2(a 



■g(a) - D^D^gia) 



n 



^ -g(b) - D m D u(S) g{b) = 0, DZ & g(0) = D+ &9 (0) and for n ± 0, g(0) = 0. From 

the boundary conditions we get at once g(a) = X^ 1 G(a) and g(b) = \~ 1 G(b). 

For n = the equation (3.12) with the boundary conditions D^y^Dy^g(a) = 
D u(y) D v(y)9(b) = and the gluing condition D~^g(0) = D±^g(0) is just the ordinary 
differential equation for a one-dimensional diffusion process that has been considered 
infinitely many times, and it has a solution for every G G C[a, b\. Let us go to the case 
ti/0. We are going to consider the intervals [a, 0) and (0, b] separately; what follows is 
about the interval (0,6]. 

Similarly to how it is done in, e.g.jl], we can prove that there exist two non-negative 
solutions and £2(2/) of the equation 

(A + ~—My) - D m D m Uv) = , < y < b , (3.13) 
A 2 (y) 

the first one increasing and the second one decreasing, £i(0) = £2(6) = 0, £1(6) < 00, 
£2(0+) = 00. The derivatives D^r^\^i(y) are increasing, -Du(£)£i(0) = 0, D^iyj&ib) < 0. 
It is easily checked that the Wronskian 



W(y) = det 



6 (y) / 

(both summands D^^^i(y) -£2(2/) and — Dq^^iv) - £i(y) are positive) does not depend 
on y: PF(y) = W > 0. 

Now we define, for y € [0, b], 



g(y) 



IT" 



£2(w) / y £i(^)-G(z)^(z) + 6(y) f &{z) ■ G{z)dv{z) 




(3.14) 



It is easily checked that Xg(y) — Ag(y) = G(y) for < y < b. 
Of course 



\m\ < 



\G\\ 

w 



£>(//> + &G7) f Uz)cB{z) 

./y 



Let us check that this goes to as y — > 0+. 
We have: 



\ + n 2 /X 2 (z) 



(3.15) 
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so the first summand in the brackets in (3.15) is less or equal 



c ,~ Du{y)tm-Du{y)Zm Uy)-Du(y)Uy) ^ w 

&KV) ; — 7F7T, - = — ; — ~ < 



minjA + n 2 /X 2 {z)\ minjA + n 2 /X 2 (z)] minjA + n 2 /X 2 {z)\ 

0<z<y ' 0<z<y 0<z<j/ 

and it goes to zero as y — > 0+. 

The second summand in (3.15) is less or equal 



niin [A + n 2 /X 2 (z)] min [A + n 2 /A 2 (z)] 

y<z<c c<z<b 

where y < c < b. The first term in (3.16) is less or equal 

-&(y) ■ D*($)&(y) < w 



jnin [A + n 2 /X 2 (z)] ~~ jnin [A + n 2 /X 2 (z)] 

y<z<c y<z<c 

and it can be made arbitrarily small by choosing a positive c close enough to 0. The 

second term in (3.16), for a fixed c > 0, converges to as y — > 0+. So we get that 

Jim g(y) = 0. 
j/->-0+ 

Now we are going to find D„^^(0+). We have: 



Du(y)g{y) 



W 



Du (F) Uv) f Uz)-G{z)dv{z) + D^y ) Uv) [\i(z)-G(z)du(z) 

Jy JO 

(3-17) 

The first integral here is equal to / + / , and it is not greater than 

Jy Jc 

\\G\\-[t 2 (y)-v(c)+Mc)-v(b)] , 

and the first summand is not greater than 

\\G\\/W ■ [W ■ v(c) + 6(c) • v(b) ■ D^Uv)] ■ 

By choosing c € (0, b) close enough to we make v(c) arbitrarily small; and we 
know Du(y\£i(y) — >■ as y — >■ 0+. So the first summand in (3.17) goes to as y —¥ 0+. 
The second summand in (3.17) does not exceed in absolute value 

IIGII • ■ \D^ 2 (V)\ ■ v(y) < \\G\\ -W-v(y)^0 (y 0+) . 

Now we are looking for the solution g(jj) of the equation (3.12) with the boundary 
conditions under this formula in the form g{y) = g(y) + C -£i(y). For the undetermined 
coefficient C we get one linear equation, and it does have a solution since / 0. 
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The same way we get, for n ^ 0, a solution g(y) for y < with g(0—) = 
Du(y)g(0-) = 0, g(a) = ^(a). 

So we get a solution / € -D(A) of the equation (3.11) for every function F(8,y) = 

N 

e in6 -G n (y), G n (y) € C[a, b], such that G„(0) = for n ^ (we take /(o) = G (0)). 

7l=-JV 



The set of such functions is dense in C(£) so that the closure operator (XI — A\ D ^) 1 
is defined on the whole C(£) which finishes the proof. □ 



Let q t be the Markov process corresponding to A\ D ^, whose existence was proved 
in Lemma 3.1. We prove the following 

Theorem 3.1. As e 1 0, for fixed T > 0, the process qf = ir(qf) converges weakly 
in the space C[ 0i t](<£) to the process q t . 

The proof is again based on an application of Lemma 2.2. 

Proof of Theorem 3.1. Making use of Lemma 2.2, we take the metric space 
M = S 1 x [a — 1, 6+1] with standard metric. The mapping Y = 7r. The space Y(M) = £ 
is endowed with the metric d, defined as follows. For any two points yi) and (82, 2/2) 
on £ with yi, ^2 having the same sign we let d((9i,yi), {82,1)2)) be the Euclidean distance 
between points (\yi \ cos^i, \y\ \ sin^i) and (\y2 \ cos^2 ; I2/2I sin ^2) in if Vi an d 2/2 have 
different sign we take d((9\, y±), (92, 2/2)) = d((^i,yi),o) + d(o, (^2, ^2))- With respect 
to this metric the space £ is a complete separable metric space. We take the process 
(Xf,~P £ x ) as q\ and the process (yt,P y ) is taken as q t . 

For the uniqueness of solution of martingale problem we set the space ^ be the 

N 

space of all continuous functions on £ which has the form F(8,y) = ^ e md • G n (y), 

n=-N 

G n € C[a, b] is continuously differentiable inside [a, 0) and (0,b], also G n (0) = for 
We take f(o) = Gq(0). It is proved in the proof of Lemma 3.1 that the equation 
Xf — Af = F always has a solution / € D C D(A) for all F € \£ and A > 0. The space 
.D contains those functions / € C(£) that are bounded and are three times continuously 
differentiable inside £+ = {(8, y) £ £ : a < y < 0} and t~ = {(8, y) € £ : < y < b}. 

We will state pre-compactness of family of distributions of processes gf in Lemma 
3.2. What remains to do is to check that for every compact K C £ and for every / G D 
and every A > we have 



E 9o 



00 

e- Xt [Xf(7r(q!)) - Af(7r(ql))]dt - f(n(q )) 









as e I uniformly in q € tv 1 (K). The proof of this is essentially the same as the 
proof we did in Lemma 2.5, based on the following auxiliary Lemmas 3.9 (for the proof 
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of convergence for processes near o) and 3.10 (for the proof of convergence for processes 
away from o) and the auxiliary Lemmas 2.9 and 2.10 (for the estimates on the exit times, 
notice that the stopping times a n and r ra we will work with in this section are essentially 
the same stopping times that we worked with in Section 2 since we are discussing about 
a model problem). We omit the details in the proof. □ 

Let k be a real number with small absolute value. Let G(k) = {(9,y) G S 1 x 
[a - 1, b + 1] : a - 1 < y < -1 - k or 1 + n < y < b + 1}. Let C + (k) = {(6, y) G 
S 1 x [a - 1, b + 1] : y = 1 + k} and C~(k) = {(9, y) G S 1 x [a - 1, b + 1] : y = -1 - «}. 
Let C(/e) = C + (k) U C _ (k). Let <5 > <5' > be small. We shall introduce a sequence of 
stopping times ro < oo < t\ < o\ < T2 < 02 < ••• by 

r = , cr n = min{t > r n , G , r n = min{t > cr n _i, q\ £ C{5')} . 

This is well-defined up to some (k > 0) such that 

{yt+a k hits a — 1 or b + 1 before it hits — 1 — 6' or 1 + 8') = 1 . 

We will then define r^+i = min{£ > : yf = a — 1 or 6+1}. And we define 
Tk+i < crk+1 = Tk+i + 1 < r k+ 2 = T fc+ i + 2 < c7 fc+2 = T fe+ i + 3 < ... and so on. 

We have lim r n = lim a n = oo. And we have obvious relations q £ T G C(5'), 
q e (Tn £ C(<5) for 1 < n < k (as long as k > 1, if k = the process may start from G(<5) 
and goes directly to S" 1 x {a — 1} or S 11 x {6+ 1} without touching C{5') and is stopped 
there, or it may start from S 1 x (—1 — 5,1 + 5), reaches C(5) first and then goes directly 
to S 1 x {a — 1} or S 1 x {b+ 1} without touching C(5') and is stopped there). Also, for 
n > k + 1 we have = q% n G S 4 x {a - 1} or S 4 x {6 + 1}. If q £ = q G G(<5), then we 
have (To = and t\ is the first time at which the process q\ reaches C(5') or S 1 x {a — 1} 
or S 1 x {6+1}. 

Note that these stopping times are the same as those defined in Section 2 since our 
process yf is essentially the process q £ in Section 2. 

The pre-compactness of the family {qf} £ >o in C[ 0i t](£) for < T < oo is proved 
in the same way as in the one-dimensional case. We shall make use of the technical 
Lemma 2.3 with ql' S and ql replaced by q^' S and ql and the space C[o,t](^) instead of 
C[ 0j t]([«; b]). We omit the proof of the next lemma. 

Lemma 3.2. The family of distributions of {qf } £ >o is pre-compact in C[ T ](£). 
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The next few lemmas establish the estimates on the asymptotic joint law of the 
processes (yf,9f) at first exit from a small neighborhood of the domain within which 
the friction vanishes. This is the key part to the proof of Theorem 3.1. 

Let 8" > be small. We consider the process q\ starting from q\ = q G S 1 x 
[— 1 — 8', 1 + 8']. Let us introduce another sequence of stopping times a± < j3\ < a 2 < 
(3 2 < ... < a n{e) by 

cki = min{0 < t < a : q\ G C(0)} , fii = min{cti < t < a : q\ € C(-S")} , 

and for k > 2 we define 

a k = minl/Jfc^i < t < a : q\ € C(0)} , (3 k = min{a k < t < a : q\ G C{-5")} . 

Here we take the convention that the minimum over an empty set is oo. The 
number n(e) is a non-negative integer-valued random variable such that a n u\ < oo and 
Pn(s) = oo. If a± = oo we set n(e) = 0. 

Lemma 3.3. For q G G(5') we have 




(3.18) 



Proof. If 1 < y £ = y < 1 + 5' 



we have 



Pq (ai < oo) 



u £ (l + 6) -u £ {y) > 
u e (l + 5) -n £ (l) - 



u £ (l + 5) -u e {l + 8') 
u e {\ + 5) -u e {l) 



= 1 - 



u{5') + sS' 
u(8) + e5 



If — 1 — 8' < yl) = yo < — 1 we have 




u e (y)-u e (-l-8') u £ {-l-8')-u £ {-\-8) 
~ u e {-l)-u e {-l-8) ~ u e {-\) -u £ (-l -8) 



= 1- 



+ e5' 
-u(-8) + e<5 



If — 1 < ?/o = Vo < 1 we have Pq () (cti < oo) = 1. □ 



Lemma 3.4. For q G G(5') we /iaue 



Pq o 0#i < oojcti < oo) > 1 — max 



5(<5) + e(<5 + 8") ' - 



e<5' 



u{-8) + e{8 + 5") 



) 



(3.19) 



Proof. If y s a = 1 we have 
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Pq (/3i < oo|ai < oo) 
If y e ai = — 1 we have 



u £ {l + 8) -u £ (l) e5" 

u £ (l + 8) -u £ (l-8") ~ ~ u{8) + e(6 + 8") 



PqJfii < oo|ai < oo) 



u £ (-l) -u £ {-\ -8) 



= 1 - 



eS" 



u £ {-l + 8") -u £ {-l -8) ' -u(-S) + e(8 + 8") 



□ 



Let M(e) — > oo as e 1 be an integer. The exact asymptotics of M(e) will be 
specified later. We prove 

Lemma 3.5. For q G we /tawe 



P qo (n(e) > M(e)\ai < oo) > 



1 — max 



e5" 



sd" 



u(S) + e{5 + 5") ' -«(-<$) + e(S + <5") 



1 M(e)-1 



(3.20) 



Proof. This is because trajectories of q\ between times ai < t < oti+i are inde- 
pendent and by iteratively using Lemma 3.4 we get the desired result. □ 



Lemma 3.6. We have 



. ' 5" x 5 
oti+i - & > e 2 ( — 



(3.21) 



with Hi being i.i.d. positive random variables with E(i7«) 4 < oo fori = 1,2, ...,n(e) — 1. 



y! 



Proof. This is a result of the Holder continuity of the standard Wiener trajectory 
\Wt — W s \ < Hi\t — sj 1 / 5 and the fact that between times < t < a.i+\ the process 

is a time-changed Wiener process -Wt traveling at least a distance of 5" . □ 

e 

Let us define an auxiliary function 



n(e,8, 8',8",M{e)) 



= 2 



1 



1 — max 
2 max 



s5" 



e8" 



M(e)-1 



u{8) + e{8 + 8") ' -u{-8) + e(8 + 5") J_ 
u{8') + e8' -u(-S') + £8' 
u(8) + £8 ' -u(-S) + e<5 



+ 
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Lemma 3.7. For q € G(S') and for some A > 0, k > and C > 0, there exists 
£q > swc/i that for all < e < Eq, for any < 9\ < 82 < 27r roe /iave 

P qo G [*i , 2 ] , ^ = 1 + 5) - ^T^o = 1 + <*) 
< Cexp(-A(5") 5 KM(£)) + 20(e,<5, <$',<$", M(e)) 

anc? 

P 9o (^ e [0i,0 2 ],<4 = -1 - 8) - ^1p, o (^ = -1 - 5) 
< C exp(-^(<5 / ') 5 KM(e)) + 2«(e, 5, <5', 5", M(e)) . 

Proof. As we have 

we set T £ (t) = / — — — -77. Using Lemma 3.6 for g n £ G(S') the random time 

T £ ((To) can be estimated from below by 



r° d, i i "l^ 1 1 



(If n(e) = 0,1 the sum is supposed to be 0.) 

And we also notice that the random time T £ {ao) only depends on the behavior of 
the process y\ and is therefore independent of the Wiener process W/ in the stochastic 
differential equation xf = — — W} (see (3.2)). For the same reason the random 

variables y^ Q , n(e) and a\ are of course also independent of W} . 

As we have the elementary inequality (e 1 ) (E(i^) 4 ) 1 / 4 > (EiJj) > 

1, we have, by Strong Law of Large Numbers 

M(e)-1 



for some constant c > 0. (We can always assume that Hi is uniformly bounded from 
below by a positive constant so that JJjy^, 1 ^ 00 an( ^ we can a PPly SLLN.) 

Now we see that we can find some £0 > such that for all < £ < £q we will have 

Pq (r £ ((7o) > (5") 5 KM(£)\n(£) > M Or) , Ql < 00) = 1 

for some constant k > 0. 
This gives 
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P 9o (T £ (a ) > (S") 5 KM(e),y £ ao = 1 + S\n(e) > M(e), ai < oo) 
= P 9() (^ = 1 + S\n(e) > M(e), ai < oo) . 

Recall that we have 9% a = x% Q mod 2ir = W^, e ,^ mod 2ir. Using this, the inde- 
pendence of T e (o"o), y^-o' a i an d n ( e ) w ith W/, and the above estimates we have, as 
< e < eo, that 



P 9o (0L, € [0i,0 2 ],z4 = ! + %(£) > M(e),ai < oo) 



/•oo 

y P qo (r £ ( CTo ) Gdt,^ = l + %(e) >M(e),ai < oo)P q() (^ 1 mod 2tt G [0i , 2 ] ) 

P qo (T £ (a ) e dt,y e ao = 1 + S\n(e) > M(e),ai < o^P,.^ 1 mod 2^ G [<9i,^ 2 ]) 



'(5") 5 AM( £ ) 

Since we have the exponential decay 

P(W t 1 mod27r € [0i,0 2 ])- 
for some C > and A > 0, we could estimate 



2vr 



< Cexp(-Ai) 



!P 9o (^o G [^'^] , <4 = ! + %(£) > M(e), ai < oo)- 



6*2 — 9\ 



2tt 



P «?o(^o = 1 +< 5 K £ ) ^ < °°)| 



< Cexp(-A(5"fKM(e)) 

for < e < Eq. 

Notice that we have, by using Lemmas 3.5 and 3.3, 



\P q Mo e [OiM,y% = 1 + &) ~P q Mo G , ^ = 1 + S\n(e) > M(e), ai < oo)| 

= |P 9o (^ G [OiMtio = ! + %(£) > M(e),ai < oo)P(n(e) > M(e),ai < oo) - 

P,o(^o € , ^ = ! + %(£) > M( £ ), ai < oo)| +P 9o (n(e) < M(e)) + P 9o (ai = 



< 2(P 9o (n(e) < M(e)) + P, (ai = oo)) 

< 2(P 9o (n(e) < M(e)|ai < oo) + 2P Qo ( ai = oo)) 

e5" e5" 



< 2 



1 



1 — max 
2 max 



-i M(s)-1 



u{5) + e{5 + 5") ' -u{-5) + e((5 + <5") J _ 
g(<5') + £<5' -u{-5') + 
u(<5) + e5 ' -u(-<5) + £(5 



+ 



Q(e,5,5',5",M) . 
By the same argument we can estimate 



2vr 



2vr 



< Sl(e,6,6',8",M) 
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Summing up these estimates we have 



< |P 9o (^ G [6>i,<? 2 ], y £ ao = 1 + 6)- P 9o (^ G [0i,0 2 ] , 24 = 1 + S\n(e) > M(s), a, < oc)\ + 
|P,o(^ G Pi.fc] , Va = 1 + S\n(e) > M(e), ai < oo)- 

P qo (f ao = 1 + S\n(e) > M(e) iai < oo)| + 



< 2n(e, 5, 5', 5", M) + Cexp(-A{5") 5 KM(e)) , 
as desired. The other inequality is established in a similar way. □ 



Combining Lemma 3.7 and Lemma 2.7 we can have 

Lemma 3.8. For q G G(5') and for some A > 0, k > and C±,C2 > 0, there 
exists Eq > sitc/i i/iai for all < e < £q, for any < 9\ < 62 < 27T we /iawe 



p, (^,G[<9i,e 2 UL = i + «y) 



#2 - 0i 5(0) - 5(-5) 



2vr 2(5) - 5(-5) 



and 



< d eM-A(5") 5 *M(e)) + 2«( e> 5, 5', 5", M(e)) + - 

9 2 -0i 5(5) -5(0) 



(3.22) 



£ 

&0 



-1-5) 



2vr 5(5) - 5(-5) 

< C 1& M-A{5"fKM{e)) + 2n(e,S,S',S",M(s)) + ^^||±p ee p(e ) 



(3.23) 



Now let us specify the asymptotic order of M(e) — > oo, 5 = 5(e) — >• 0, 5' = 5'(e) — > 
and 5" = 5"(e) — >• as £ J, 0. Since for < k < 1 we have the elementary estimate 
1 — (1 — tx) n = k(1 + (1 — k) + ••• + (1 — k) ti_1 ) < Kn we can estimate 



ft(e,5, 5',5",M(e)) 



< 2 



M(e) • max 



eS" 



e5" 



5(5) + e(S + 8") ' -5(-5) + e(S + 5") 
5(5') + eS' -u(-8') + s5' \ 

- niax| ~(5)+ £ 5' -5(-5) + £ 5 J • 



+ 



We shall choose 5" = 5"(e) << 5 and M(e) such that the requirements of Lemmas 
2.6, 2.7 and 2.8 hold. At the same time, we need 

1 



(5"fM(e) > In 



(5(5) - u(-6))< 



(3.24) 
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and 



M(e) 



(3.25) 



To this end we let M(e) = In ( ^ ) and 5" = 



1\ .„ / (u(S)-u(-S)^ U (u(S)-u(-S)^ 2 ^ 



hi: 



. At 



the same time we keep our asymptotic order of choice of e, 5 and 5' as in Section 2. This 
means that we need 

4/5 1 / 1 \ 2 

<(5(*)-5(-«)) 2 . 



e In 



In 



v eyy u(5) — «(— 5) \u(S)—u(—S) 

It could be checked that this is possible to make (3.24) and (3.25) to hold. We formulate 
this as a corollary. 



Corollary 3.1. Let q Q G G(S'). Under the above specified asymptotic order we 
have, there exist eq > such that for any < e < Eq we have 

< 2 - 0i 5(0) - 5(-5) 



P qo (^ o e[0 1 ,0 2 ], 2 4 = l + <5)- 



2vr 5(5) - u{-8) 



< C ■ (u(5) - u(-5)Y , (3.26) 



P, (^€[01,02], = -!-<*) - 



5(5) - 5(0) 



2vr 5(5) - 5(-5) 



< C-(5(5)-5(-5)) 2 . (3.27) 



Lemma 3.9. For any q £ G(5') and for any p > there exist £q = £o(/°) such that 
for any < e < Eq, for any f € D(A) we have, for some K > 



E q /(7r(<4)) - /( W (g))| < K(5(5) - 5(-5)) 2 



(3.28) 



Proof. We have, using Corollary 3.1, that 



V q f{ir{q% Q ))-f{ir{q))\ 
E,/(^,7r(^))-/(7r(g))| 

J f(9,6)P q (e* ao ede,tf ao = l + 5) + J /(^-«5)P q (0- o Gd0,^ o = -l-«5)-/(7r(Q)) 

J_ 5(0) -5(-5) ^ Ja t 1 /- 27r 5(5) -5(0) 
27 io 5(5) - 5(-5) 



< 



5(0) - 5(-5) 



2tt Jq u(S) — u(—5) 



f(9,-S)dB-f(ir(q)) 



+ K 1 {u{5)-u{-5)f 



(f(0,6)-f(o))dB 



< 



2tt Jq u{8) — u(—6) 

h r w^wT) im - m - 5))d9 + (/(0) - f ^ q))) 

(5(0) - 5(-5))(5(5) - 5(0)) / 1 [*« f(9, 5) - f(o 



5(5) - 5(-5) \2tt J 5(5) - 5(0) 

|/(o) - f(*(q))\ + ^(5(5) -5(-5)) 2 
< K(u(S) - 5(-5)) 2 



-de 



i 



2tt 



+ Ki(5(5) -5(-5)) 2 
/(o)-/(0,-5) 



2vr 7 5(0) - 5(-5) 



-de 



+ 
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for some K\ > and K > 0. We have used the gluing condition (3.10) and our specified 
choice of asymptotic order of S, 5' and e. □ 



Lemma 3.10. We have, as e,5,5' are small, for q G G(5), that 



E 



f 1 e - At [A/(7r(q?)) _ A/(*r(g?))]dt + e^/W^)) 



/«<*))) 



< (S(5)-S(-«)) 2 . 
(3.29) 

The proof of this Lemma is essentially the same proof in Lemma 2.6 modified into 
a two-dimensional version and we omit it. 



Finally we would like to mention that our boundary condition given in this section 
also appears naturally in other model problems. As an example let consider the following 
system: 

ft i 
x\ = / — — dW} , 

1 J X(y!)+e 1 (3.30) 

Vt = \w t 2 \ • 

Here A(«) is a smooth function on R + that vanishes at and is strictly positive 
in (0, oo); W} and are two independent standard Wiener processes on R. Let 
the process z\ = (xf,yf) on R x M + be stopped once it hits the boundary {(x,y) G 
R 2 : y = R} for some R > 0. Let 9f = x% mod 2vr. Let tt : S* 1 x R + -»• M 2 be 
the mapping defined by iv(6,y) = (y cos 6 , y sin 8) . For each fixed e > 0, the process 
iuf = (0f, yf) is a diffusion process on S" 1 x [0, R] with normal reflection at the boundary 
{(0,y) : y = 0} and is stopped once it hits the other boundary {(0,y) : y = R}. Let 
mf = Tt{wt) (i.e., we glue all points {(0,y) : y = 0}). The process mf moves within 
the disk B(R) = {m £ R 2 : |m| K 2 < R} and is stopped once it hits the boundary. In 
general, this process is not a Markov process. But we expect that, as e \, 0, this process 
will converge weakly to a Markov process wt on B(R) with generator A and the 
domain of definition D(A), defined as follows: The operator A at points (9,r) (we use 
polar coordinates, that is, a point (x,y) € R 2 is represented by (r cos 9,r sin 9)) with 
r 7^ is defined by 

The domain of definition D(A) of the operator ^4 consists of those continuous functions 
/ on B(R) for which Af(9,r) is defined and continuous for r / 0, the derivative in r 
being continuous; such that finite limit 

lim %(9',r) (3.32) 
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exists; 



exists and does not depend on 8; 



and 

I-2-k 







lim Af(9',r) (3.33) 



lim A/(0',r) = O; (3.34) 



lim ^-(9',r)d0 = . (3.35) 
»0,r->o+ ar 



We define, for / G A/(0, fl) as the limit (3.34) and Af(0) as the limit (3.33). 

The weak convergence of wf to wt in C[Q t T](B(R)) described above shall be a result 
of fast motion x\ running at the local time of the slow motion yf on the boundary 
{(x, y) G R x R + : y = 0}. The proof of this result shall follow the same method of this 
section. 



4 A conjecture in the general multidimensional case 

In this section we give a conjecture in the general multidimensional case. Consider 
the general multidimensional problem (1.8), and for brevity assume that &(•) = 0. That 
is, the system has the form 

Let us work in a large closed ball B(R) = {q G M. d : \q\^d < R} for some R > 0, 
i.e., the process q\ is stopped once it hits dB(R). Suppose the friction A(«) is smooth 
and X(q) = for q in some region G C B(R) while X(q) > for q G B(R)\[G] (here [G] 
is the closure of G with respect to the Euclidean metric in M d ). The domain G C B(R) 
is assumed to be simply connected and to have a smooth boundary dG. 

Let £ be a topological space consisting of all points in B(R)\[G] and one additional 
point o. The topology of £ contains all the open subsets (in standard Euclidean metric) 
in the induced topology of B(R)\[G] and all the open neighborhoods of [G] in B(R) as 
the open subsets of £ containing o. Let us consider a projection 7r : B{R) — > <£ defined 
as follows: for points q G B(R)\[G] we have 7v(q) = q and for points q G [G] we have 
7t(q) = o. Under the above defined topology for £ the mapping 7v is continuous. Let 
qf = 7r(gf) be a stochastic process with continuous trajectories on <£. 

Our conjecture is about the weak convergence, as e 1 0, of gf to some Markov 
process q t on C Below we give our definition of the latter process but we point out that 
we are not clear about the existence of it. Our generator and boundary condition for this 
process is more or less in the spirit of martingale problems (see, for example, Ch.4]). 
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To ensure the uniqueness of solution of martingale problems we need the existence of 
solution in a nice space of the corresponding PDE with the specified boundary condition. 
We are not clear about this yet. 

The process q t is identified by its generator A with domain of definition D(A). For 
a function f(q) on £ that is continuous on <£ and smooth for q ^ o, \q\^d < R we define 

= 2A3(g) + 2A^) A/( ^ ' (42) 

and at the points o and q with |q|]gd we define the values of Af as the limits of the 
values given by (4.2) (assuming these limits exist). The domain D(A) is defined as the 
set of functions / such that Af(q) = for \q\^d = R, the generalized normal derivative 

Z3 g /(g)=lim /(g + gw y ) -/ (0) (4.3) 

exists for all q € dG, where n(q) is the vector of the outside normal to dG, and u(q) is 
some function defined in a neighborhood of dG with lim u{q) = 0; and 

7r(q)^0 

/ D u f{q)da{q) = . (4.4) 
JdG 

Here da(q) denotes integration with respect to the surface area on dG. 

Conjecture. The process (f t = 7r(qf) converges weakly in the space C[q ) t](£) as 
e 10 to the process q t described above. 

A further conjecture: we can define the function u as 

u(q + 8n{q))= I \{q + sn{q))ds (4.5) 
J o 

for q € dG and S > sufficiently small. 
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